Singularities in Speckled Speckle 
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Speckle patterns produced by random optical fields with two (or more) widely 
different correlation lengths exhibit speckle spots that are themselves highly 
speckled. Using computer simulations and analytic theory we present results 
for the point singularities of speckled speckle fields: optical vortices in scalar 
(one polarization component) fields; C points in vector (two polarization 
component) fields. In single correlation length fields both types of singularities 
tend to be more-or-less uniformly distributed. In contrast, the singularity 
structure of speckled speckle is anomalous: for some sets of source parameters 
vortices and C points tend to form widely separated giant clusters, for other 
parameter sets these singularities tend to form chains that surround large 
empty regions. The critical point statistics of speckled speckle is also anoma- 
lous. In scalar (vector) single correlation length fields phase (azimuthal) 
extrema are always outnumbered by vortices (C points). In contrast, in 
speckled speckle fields, phase extrema can outnumber vortices, and azimuthal 
extrema can outnumber C points, by factors that can easily exceed 10 4 for 
experimentally realistic source parameters. 
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Studies of the singularities of random optical fields 
(speckle patterns) normally involve fields with a single 
correlation length A[Diai31lllSllSllZllSllSlin3iniIiaiI31 
HI US HS1 HZ1 HH1 US EH- Here we study random fields 
with two widely different A, finding anomalous spatial 
distributions for the singularities, and anomalous statis- 
tics for them and their associated critical points. These 
interesting anomalies can be controlled by varying source 
parameters, and although studied here by computer sim- 
ulation and analytic theory, they should also be amenable 
to experimental observation. 

The presence of two distinct A in a random optical field 
is easily determined visually, because the major speckle 
spots in the speckle pattern are themselves highly speck- 
led, Figs, [I] and [2] We call such patterns "speckled 
speckle" [21]. 

Speckled speckle, which arises from source distribu- 
tions with widely different length scales, is often the 
norm. It appears in an uncontrolled form in many scat- 
tering experiments due to parasitic scattering, but be- 
cause the detector resolution is often insufficient to re- 
solve the small speckle spots, the presence of speckled 
speckle may not be noticed. But it makes its presence 
felt in discrepancies between experiment and theory, and 
is, for example, the most common cause of the failure 
to obtain negative exponential statistics for the inten- 
sity. Due to interference, parasitic fields with intensities 
less than 10 -8 of the main beam can produce measurable 
effects. 

Here we consider controlled situations in which the 
long correlation length A a and the short correlation 
length, A5 and their associated intensities I a and I bl are 



chosen to have values that are experimentally realistic. 

Commonly used rotationally symmetric (isotropic) 
source distributions S(r'), and their associated auto- 
correlation functions W(r), for circular Gaussian fields 
PQ with a single correlation length are: (i) a Gaussian, 
S G (r f ) = exp[-r /2 /(2p) 2 ], W G (r) = exp(-/€ 2 p 2 r 2 ), (m) 
a disk, Sn(r f ) = 0(r' — p), Wo(r) = 2J 1 (npr) /(npr), 
and (in) a thin ring, modeled as SR(r f ) = S(r f — p), 
Wr(t) = Jo(ftpr). Here J n is a Bessel function of integer 
order n, and n = 2ir/(\Z), where A is the wavelength 
of light, and Z is the distance to the screen in the far 
field on which the speckle pattern is measured [T]. In 
what follows we set for convenience k = 1. In each of the 
above three cases A ~ 1/p. 

Controlled fields with two different correlation lengths 
and known correlation functions are easily generated us- 
ing linear combinations of two of the above source func- 
tions with different values for the parameter p. There 
are nine such combinations; here we consider the three 
in which both source functions are the same - two Gaus- 
sians, two disks, and two rings. Writing: (i) for the Gaus- 
sians, S a (r') = exp[-r /2 /(2a) 2 ], W a (r) = exp(— a 2 r 2 ), 
S b (r') = exp[-r /2 /(26) 2 ], W b (r) = exp(-6 2 r 2 ); (ii) for 
the disks, S a (r') = B(r' - a), W a (r) = 2J 1 (ar)/(ar), 
S b (r') = S(r' - 6), W b (r) = 2J 1 (br)/(br)- and (in) 
for the rings, S a (r') — 5(r f — a), W a (r) — Jo(ar), 
S b (r f ) = S(r' — 6), Wb(r) — Jo(br); we have for our 
compound source function S ab (r f ) and its scalar far-field 
autocorrelation function W a b(r), 

S ab (r f ) = I a S a (r') + hS h (r') (la) 
W ab (r) = [W a (r)+ KW b (r)]/(l + K) , (lb) 



where for the Gaussians and the disks K = b 2 Ib/ '(a 2 I a ), 
and for rings K = blb/(al a ). Here A a ~ 1/a, A5 ~ 1/6. 

The point singularities of generic scalar (one polar- 
ization component) fields are optical vortices [2j [22] , 
whereas for generic vector (two polarization component) 
fields the point singularities are C points - isolated points 
of circular polarization embedded in a field of ellipti- 
cal polarization [12] [T3l [22]. It is convenient to treat 
C points as vortices of the complex Stokes field 612 = 
Si + iS2 [23] [24], where in terms of Cartesian field com- 
ponents E x and E y , and right (R) and left (L) circularly 
polarized components, En and El, 

51 = \E X \ 2 - \E y \ 2 = 2Re(E* R E L ), (2a) 

5 2 = 2Re(E*E y ) = 2Im(E* R E L ). (2b) 

Here Cartesian and circular field components are related 
by 

Er = (E x - iE y )/2^ 2 , E L = (E x + iE y )/2^ 2 (3) 

C points are also optical vortices of Er and of El, be- 
cause at a point where one of these fields vanishes, as 
happens at an optical vortex, the optical polarization is 
circular, whereas elsewhere in the field when both compo- 
nents are nonzero the polarization is in general elliptical. 
Thus the C point density is twice the vortex density. Fi- 
nally, for vortices, scalar fields a and b have the same 
polarization, whereas for C points in vector fields a and 
b are orthogonally polarized. 

These essential preliminaries completed, we turn now 
to our results for speckled speckle in what may be called 
the perturbation regime, i.e. the regime in which Ib/I a <^ 
1. Results for other regimes will be reported on sepa- 
rately. 

In Fig. 1, which shows scalar fields, we compare normal 
speckle (lb = 0) with speckled speckle (lb > 0). As can 
be seen, the vortices of speckled speckle are found only 
in regions surrounding the a field vortices shown in Fig. 
[TJb). The reason for this is that at a vortex the total, 
i.e. composite, field amplitude must vanish. But because 
lb <C I a , the destructive interference between the a and 
b fields required for vortex generation can occur only in 
regions where I a = 7^, i.e. where I a is very small. This 
leads to isolated giant clusters of speckled speckle vortices 
near the positions of the a field vortices. Such a cluster, 
containing 29 positive and 28 negative vortices is shown 
in Fig. [lje). As shown in Figs. [TJf) and[TJg), vortices 
are absent in other regions. 

An anomalous feature of speckled speckle phase maps 
is a granularity not seen in normal speckle phase maps; 
not only the intensity, but also the phase of speckled 
speckle is speckled. This granularity, which is present 
throughout the phase map, and which can be seen most 
easily in Fig. [ljg) , corresponds to large numbers of phase 
maxima and minima, i.e. phase extrema. 

In normal speckle phase maps vortices V always out- 
number extrema E [6] [14]. For a single Gaussian E/V = 
0.414, for a single disk E/V = 0.0887, and for a sin- 
gle ring E/V = [14 . E/V can be expressed in terms 
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Fig. 1. Scalar speckle. Shown are simulations for disks 
with a — 1, b — 100. In (a) and (b) lb = 0, which yields 
familiar single A speckle with A ~ 1. In (c) - (g)Ib/I a — 
2 • 10" 6 (K = 0.02), which yields speckled speckle with 
A a ~ 1, Ab ~ 0.01. (a) Intensity of single A speckle, (b) 
Phase map corresponding to (a), (c) Intensity of speckled 
speckle, (d) Phase map corresponding to (c). (e) - (f) 
Forty-fold enlargements of the regions in (d) pointed to 
by the arrows labeled e, f, and g. Here and throughout 
intensity (phase) is coded increasing (0 to 2ir) black to 
white, and positive (negative) vortices are shown by filled 
white (black) circles with black (white) rims. 



of the moments M n , n = 2,4, of the source function 
S [14] . For speckled speckle we obtain for these moments 
M n = C n (a n + Kb n )/(l + K), where for two Gaussians 
(2G), C n = 2 n (n/2)!, for two disks (2D), C n = 2/(n + 2), 
and for two rings (2R), C n = 1. Writing p = a/b, when 
K = p 2 , E/V attains a maximum value that for small p 
is given by 

(E/V) max =g[(3)^ 2 p 2 ]-\ (4) 

where for 2G, g = 1/2, for 2D, g = 1/3, and for 2R, 
g = 1/4. In principle (E/V) max can be made arbitrarily 
large; for the parameters a and b in Fig. 1 we have for 
2G, 2D, and 2R, respectively, (E/V) max = 2887, 1924, 
and 1443, corresponding to enhancements over normal 
speckle by factors of ~ 7, 000, ~ 22, 000, and oo. 

C points show similar anomalies. At a C point \E X \ = 
\E y \, so also speckled speckle C points can be found only 
in dark regions of the a field, here E x . In Fig. 2 we show 



3 



"7T 



M 



? a* • ~ o 







• 














(d) 


(e) 


(f) 



Fig. 2. Vector field speckled speckle. Shown are simu- 
lations for disks with a = 1, b = 10, Ib/I a — 3 • 10 -3 
(K = 0.3). (a) Intensity, (b) Phase $12 of complex 
Stokes field S12. (c) Enlargement of the black- white band 
on the horizontal centerline in (b). (d) - (f) Enlargement 
of the region beneath the black- white band in (b). (d) 
Phase $12 of 5i 2 . (e) Phase —$r of E R . (f) Phase § L 
of E L . 



speckled speckle C points for a value for Ib/I a that is 
much larger than the corresponding value in Fig. [I] Now 
C points are no longer confined to the darkest regions 
near a field vortices, but can populate dark areas between 
a field speckle spots, typically forming chains that tend 
to surround these spots. An example of such a chain, 
which contains 28 positive and 29 negative C points, is 
shown in Fig. [2|c). 

Also speckled speckle C point phase maps contain nu- 
merous extrema that are absent in normal speckle C 
point maps: in obvious notation E/C = \E jV . This 
result follows from the fact that, as shown below, the 
density of phase extrema is essentially the same in both 
scalar and vector speckled speckle, whereas the density 
of C points is twice the vortex density [I2j [I3j [22] . 

From Eqs. (3) we have E R = E x (l — iE y /E x ) 
« E x exp(-iE y/ E x ), E L = E x (l + iE y/ E x ) « 
E x ex.p(iE y /E x ), where in the approximate equalities we 
have used the fact that E x is the strong a field, and E y is 
the weak b field, so that E y jE x is small. Writing Si 2 — 
Ai 2 exp(z$i 2 ), E x = A x exp(i(p x ), E y = A y exp(i(p y ), 
Er — Ar exp(i<I>#), and El = Al exp(z$^), and using 
Eqs. (3), we have from Eqs. (2) § 12 = $l - §r~2(§ l - 
(p x )~ — 2(&r — (p x ). But because A x = A a ~ 1/a is large, 
except close to its sparse vortices ip x varies slowly com- 



pared to $12, and for which Ai 2 ~A^^A/ / ~ 1/6. 
Thus except close to its sparse vortices tp x can be taken 
to be nearly constant, and so the phase field $12 of vec- 
tor speckled speckle, and the phase fields <&l and — &r 
of scalar speckled speckle, have very nearly the same set 
of extrema. This unusual property of speckled speckle is 
verified in Figs. 2(d)-2(f). 

In summary, we have shown that in the perturbation 
regime studied here the singularities of scalar and vector 
speckled speckle, and the critical points associated with 
these singularities, have highly interesting, anomalous 
spatial arrangements and statistics. Singularity screen- 
ing, correlations, and other properties of these unusual 
fields in different parameter regimes are also likely to ex- 
hibit unusual new phenomena of special interest. 

Email addresses: I. Freund, freund@mail.biu.ac.il; D. 
A. Kessler, kessler@dave.ph.biu.ac.il. 
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